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Abstract 

We introduce a simple continuous model for nonequilibrium surface growth. 
The dynamics of the system is defined by the KPZ equation with a Morse-like 
potential representing a short range interaction between the surface and the 
substrate. The mean field solution displays a non trivial phase diagram with 
a first order transition between a growing and a bound surface, associated 
with a region of coexisting phases, and a tricritical point where the transition 
becomes second order. Numerical simulations in 3 dimensions show quanti- 
tative agreement with mean field results, and the features of the phase space 
are preserved even in 2 dimensions. 

I. INTRODUCTION 

Recently Hinrichsen et. al. [Q have introduced a discrete model for nonequilibrium surface 
growth that displays a first order unbinding transition in a region of the parameter space 
where there is phase coexistence between a growing and a bound surface. They used a 
restricted solid on solid dynamics, and distinguished the probability of adsorption of particles 
at the substrate from that at higher layers. This in order to take into account short range 
interactions between the substrate and the surface. 

In the same spirit we study a continuum growth model for an interface interacting with 
a substrate through a short range attractive potential. The interface separates a solid from 
a vapor phase, and is driven by a difference in chemical potential due, for example, to a 
flux of incoming particles, as in Molecular Beam Epitaxy. As the driving force increases the 
interface eventually detaches from the substrate and grows indefinitely according to KPZ |2j 
dynamics. We will focus on the nature of the transition from a bound to a moving phase. 
We resort to numerical simulations and to a mean field approach first introduced in || for 
the KPZ equation. 

A similar unbinding process of a surface from a wall has been studied by Mufioz and Hwa 
M. They considered a KPZ M equation to which they added strictly attractive or repulsive 
interactions, and identified a second order unbinding transition for which they calculated 
some critical exponents. 

We will show how the short range attractive potential changes the nature of the phase 
transition, which becomes first order for an entire range of values of the parameters of the 
model. The second order transition is recovered only for larger values of the noise, or in the 
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case of a very short ranged potential, but the critical exponent is still different from that of 
the simple repulsive wall case. 

Our model is an example of continuous nonequilibrium system with noise. In the context 
of noise induced phase transition || this kind of systems has been studied by Van den Broeck 
et al. [§,0, who showed that multiplicative noise is essential for the transition to take place, 
since it can trigger instabilities in the short time dynamics of the system. On the other 
hand systems with multiplicative noise are characterized by a transition to an absorbing 
phase ||, a state from which the system cannot escape, characterized by a null value of the 
order parameter. The phase diagram of this class of systems has been extensively studied 
in recent years [p|-|TT|j. Most of these works were concerned with second order transitions, 



but an example of first order has been studied by Miiller et al. [|12| with the help of mean 
field techniques and numerical simulations. 



II. THE MODEL 

With the term surface we identify the interface between the solid and vapor phases of 
some substance. In what follows we want to give a statistical physics description of the 
interactions that may take place between the vapor particles as they are depositing on some 
kind of substrate. 

First of all a note on terminology: an interaction is called attractive for the surface when 
there is some force that repels incoming particles, so that the surface is not growing; On the 
contrary, attractive interaction between the substrate and the particles will let the surface 
grow and eventually detach from it. 

In our model the system is described by a height field h(x) defined on a rf-dimensional 
continuous substrate, which evolves according to a Langevin equation of the KPZ type. The 
latter is chosen in order to have a system out of equilibrium. 

The Langevin equation for the field h(x) is obtained in the usual way: 

d t h(t)=D[V 2 h-(Vh) 2 ]+r- S -^ + ar ] , (1) 

where D is some diffusion coefficient, r is the driving force due to the incoming particles, 
and the nonlinear term (Vh) 2 characterizes the nonequilibrium KPZ-like || dynamics, and 
allows for the growth along the direction normal to tilted areas of the surface. This term 
breaks the up-down symmetry of the system and therefore defines a preferential direction 
for growth (see e. g. fll3"| ). 

The field t]{x,t) is a white Gaussian noise with the following properties: 

(7](x,t)) = 

(r](x, t)r]{x', t')) = 5{x - x')S(t - t') . (2) 

We model the interaction of the substrate with the surface introducing a potential V[h] 
which has the shape of the Morse potential (see Fig. |1]): 

V M [h\ = A(l - e - ph e pho f , (3) 
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where A is the depth of the well, (3 the hardness of the repulsive wall that mimics the 
substrate, and ho the position of the local minimum, which can be set to without loss of 
generality. 

Note that the driving force r can be included in the potential as a linear term in the 
height field: —rh(x). The parameter r is the essential one: if it is negative the shape of V 
is a well with an exponential wall on the small-h side, whereas if it is positive the minimum 
close to the repelling wall is only a local one (Fig. [I]). Particles incoming from above in this 
case have to trespass a potential barrier (i. e. the local maximum of V) in order to arrive 
close to the substrate. This is equivalent to a repulsion for particles, or an attraction for 
the surface. Physically this term has a precise meaning, since it is due to the difference 
in chemical potential between the solid and vapor phases, and determines a drive of the 
interface toward higher or lower values, depending on its sign. 

In the dynamics of growing surfaces this term is usually eliminated via the transformation 
h — > h — rt, which correspond to viewing the system from a moving reference frame. As 
one can see the potential in Eq. ([|) is not invariant for such a transformation, nor for 
the more general h — > h + 5h. Notice that the KPZ equation is invariant for both these 
transformations. In the present case the consequence is that we can uniquely define the 
mean position of the surface with respect to the substrate, and its state of motion. It is 
therefore meaningful to study the unbinding of the surface from the substrate. 

By inserting the Morse potential one obtains: 



equation with multiplicative noise; this suggests us the order parameter to be used in the 
search for a phase transition, namely the spatial average of ip, m = (i/j): 



Notice that we have used the Ito calculus in doing the transformation, and we have been 
able to set A and D equal to one via a suitable rescaling of space and time. 

Qualitatively we can have two possible situations in the long time limit of the stationary 
solution of the equation. Either m —>■ 0, which corresponds to a diverging value of (h), 
i. e. to a moving (growing) surface; or m ^ 0, corresponding to the surface being bound to 
the substrate at some average height. The transition between the two phases is controlled 
by the parameter r. Namely, as one can see from Eq. @, for large negative values of 
r the surface is pushed down against the substrate, whereas for large positive values the 
surface is pulled away from it. It is the balance between the force r and the repulsion from 
the exponential wall e~ 2/3h that determines the equilibrium position of the surface. Notice 
finally that in the language of multiplicative noise systems the state with m = is the 
absorbing one, since the surface can never invert its average velocity and go back toward 
the substrate. 

An equation similar to (|5]) has been studied by Munoz and Hwa in Q in the context of 
multiplicative noise processes. The system had a soft lower repulsive wall represented by a 
term —pip p+1 in the Langevin equation, the hard wall limit corresponding to p — > oo: 




(5) 



d t ip = V 2- ?/> + rip — pip p+1 + iparj 



(6) 
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The authors found a second order phase transition between a bound and a moving phase. 
From scaling arguments, the exponent that controls the vanishing of the order parameter 
close to the transition is found to be (3 T = (2 — z)/(2z — 2), where z is the KPZ dynamical 
exponent. Numerical simulations confirmed this result. 



III. MEAN FIELD APPROXIMATION 

To proceed further with the analytical solution of the model we notice that dimensional 
analysis performed on Eq. [5] along the lines of shows that the upper critical dimension 
for the depinning transition is 2. This means that mean field results should be correct 
in 3 dimensions, and provide a reasonable estimate (up to logarithmic corrections) in two 
dimensions, i. e. for real surfaces. 

First of all we perform a discretization in space, by inserting a (hyper)cubic lattice, and 
define the variables ip(]t,t) — > ip(xi,t) = ipi{t). The discretized Laplacian is: 

where the sum runs over the nearest neighbor of the site i. The mean field approximation 
we are going to use consists PJITf in the following expression for the discretized Laplacian: 

\ E u>j-i>i)- w-v>* • (8) 

This is equivalent to say that the number of dimensions is very large, or that the number 
of neighbors (coordination of the lattice) is infinite, or again that the interaction is infinite 
range. 

Equation ([5]) then becomes: 

d t ^{t) = (-r + a 2 /2 - 1)^ - 2(3 [^f +1 - ^f +1 ] + W - °AVi , (9) 
from which, going back to the field h, the discretized version of Eq. (d) becomes: 

d t hi(t) = (r + 1) + 2f3e~^{e-^ - 1) - {e~ h )e hi + a Vi . (10) 



The Phase Transition 

We obtain the probability distribution of the stationary solution of Eq. fllQl) by solving the 
associated Fokker-Planck equation. The stationary solution is: 

2V e// (fe) 

P st {h)=Ne , 

where M is the normalization factor, and the effective potential is in the present case: 

V eff (h) = -(r + l)h + e- 2ph ~2e~^ h + (e- h }e h . (11) 



This expression for the probability distribution depends on the value of the order parameter, 
which must be determined self-consistently. We have to solve therefore the following mean- 
field equation: 



m = (e- h ) 



[ e- h P st (h)dh = F(m) = -M?L (12) 
J g.J-i{m) 



, \_ / r, -(j+l)h -A,e-^ h XeT^ -^S-e h / 1Q \ 

gj(m) = / a he u ' e ^ e^ 7 e ^ , (13) 



where J = — 2(r + l)/a 2 . 

Plotting m and F(m) on the same graph, it is easy to see that Eq. fll2]) can have one, two 
or three solutions in m, according to the value of r. The possible situations are represented 
in Fig. The stability of the solutions can be investigated making use of the results of 
Shiino JOA. Namely, the solution can be stable only when F(m) crosses the line m from 
above. In our case we have that for r larger than some r C2 the only solution (therefore a 
stable one) is m = 0; for r below some value which will be determined in a moment, there 
are only 2 solutions, the stable one being m/0; finally for intermediate values of r there 
are three solutions, two of which look stable. Notice that m = remains always a solution; 
it becomes stable only when r reaches the value for which F'(m = 0) = 1. This value is 
found analytically to be: r c = This number is also the lowest value of r for which three 
solutions exist. 

The situation just described represents a typical case of first order phase transition, and 
is the same as found in [JT2]. The value of the order parameter m remains different from 
zero as long as this is the only stable solution. Then for r > r c there is a coexistence region, 
where two phases are stable, and the final state of the system is determined by the initial 
conditions. Finally there is the value r = r C2 above which the order parameter can be only 
zero. If we remind the definition of m we understand that this phase transition corresponds 
to the unbinding of the surface from the substrate, (e~ h ) — ► meaning (h) — > oo. We remark 
that with the term coexistence we do not mean that both phases are present at the same 
time, rather that depending on the initial conditions both are possible stationary state of 
the system. In Fig. |3] is shown an example of this phenomenon. In the coexistence region 
an interface which is initially close to the substrate will remain bounded to it whereas if the 
interface is initially far it will move further away from the substrate with a velocity v > 0. 



The Moving Phase 

It is also possible to determine the velocity of the surface when it is away from the substrate. 
To do this we define a field <f>i(t) for the fluctuations of h around its mean position: <pi{t) = 
hi(t) — vt. Here v is the mean velocity of the surface obtained from Eq. (|T0|): 

v = d t {h i ) = r + l-{e- h )(e h ) , 

where we could neglect the exponential terms in the potential, since they are vanishingly 
small in the moving phase. The field <f>i(t) obeys in this phase to the following Langevin 
equation: 
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d t <fii = d t hi 



(r + 1 



e- h )e h ' + ara 



The stationary probability distribution from the associated Fokker-Planck equation is as 
usual (see e. g. Hl4fl): 

P st (<t>)=Afe-^^ =Afe a ^ , (14) 
with a = -^(r + 1 — v) and 7 = ^-(e - ^). The self-consistency relation for (e~ h ) is: 

from which we have that v = r — The condition v > shows that r c = er 2 /2, as 
anticipated. In the next section we will see that for large values of the noise strength a and 
of /3, the unbinding transition becomes second order. In that part of the phase diagram we 
can define in the vicinity of the transition a velocity critical exponent 9 []13| |: v ~ (r — r c ) e . 
From the previous discussion we see that 9 = 1. 



IV. NUMERICAL FINDINGS 



To study further the phase diagram of the model in mean field approximation we solve 
numerically Eq. flT2p. In particular we have determined the value r* for which its two non- 
zero solutions coalesce with the null one. Analytically this point corresponds to a negative 
second derivative of F(m) at the critical point r c . In Figs. ^ and [5] are shown two typical 
cuts of the phase diagram at fixed values of the control parameters a and (3 respectively. 
One can see that as the value of one of the parameters increases while the other remains 
fixed, the coexistence region becomes smaller, and eventually disappears for some (3* (a) or 
u*{(3). At the same time the unbinding transition becomes second order, thus defining a line 
of tricritical points in the phase diagram (a,f3) of the system. A numerical estimate of the 
exponent (3 T governing the transition (see e. g. |IT| ) reveals that, at least for large values of 
/3, we have: 



m ~ r — T c 



r — r r 



72 



(15) 



This can be checked analytically by taking the limit for /3 — > 00 at the critical point in 
Eq. fll^p . The condition to have a second order transition is that F"(m) < at r c ; in this 
case there can be no coexistence region, since r C2 = 
we have: 



With the same notation of Eq. ([12] 



yj-i 



F (m) = - 




9j-2 _ gjgjs 
9j-i gj-i 



, 9 J 9 j -2 
9j-i J 



(16) 
(17) 



Then we make the change of variable x = 2m/a 2 e h , to get: 
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9j{m) 



^y +1 I™ dxx-^e-* 



2m 

2rn 



r(-j - 1) - F* dxx~ {j+2) e~ x 
Jo 



Since we are looking for a second order transition, we need to estimate the small m limit of 
gj(m). Therefore we can expand e~ x in the integral, and solve it explicitly to any order in 
m. At the critical point r c = a 2 /2 we find that: 



F (m) = --4 r < 



r(£ + i) W 2 

which means that the transition is always second order at large (3, and that the m-expansion 
of F{m) close to the critical point can be written as: 



F(m) = (1 + 5r)m + Cm' 



r+l 



where 5r = r c — r, and C is some constant. Solving the mean field equation F(m) — m = 
gives: 

m ~ (Sr) 2 , (18) 

in agreement with the numerical value. 

We note that the latter calculation can be repeated for the model of || Eq. (|5]), to give 
the same result, thus implying that in the mean field limit of a system with a hard wall the 
phase transition is governed by a continuous exponent. To our knowledge this finding was 
not previously reported. A different value for the j3 r exponent, namely (3 r = 1/p, is given 
p| , where the mean field approximation Eq. (§) is also compared with different results 



m 



from field theoretical calculations (j3 r = 1). However it is not clear whether the hard wall 
limit can be obtained from the aforementioned calculations. 

Another feature of the phase diagram is that in the limit of small a r C 2 -> (3/2 , but we 
were unable to derive this result analytically from Eq. [L2[ 

We check if there is any functional relation between (3* and a*. Figure |6] shows a log- log 
plot of the tricritical points together with the linear fitting. The value of the slope is in this 
case about —4, which means that (3* ~ cr* -4 , or that (3* ~ r*~ 2 . Notice however that this 
linear fitting is a fairly poor approximation, and there is no reason why we should expect a 
priori any simple functional relation between these two quantities. 

We have also checked with simulations on a finite number of sites that these results are 
the same as those obtained from Eq. fll0|) . With the help of these simulations we could 
see that hysteresis associated with the first order transition and the coexistence region does 
indeed occur. 



V. RESULTS IN FINITE DIMENSION 

To perform numerical simulation of Eq. (Q) in finite dimensions we introduced as usual 
a (hyper) cubic lattice on the space of coordinates, and used for the discretized Laplacian 



7 



the expression in Eq. (|7|). Going back to the h field we obtained the following Langevin 
equation on a lattice: 



d t h i = r + l + 2(3e-^ H {e~ phi - 1) - e h *- + a ^ ■ ( 19 ) 

To have a better accuracy we chose to apply the Heun method for solving stochastic differ- 



ential equations (see e. g. |lq , |l7| ). We run the simulations long enough for the system to 
reach its stationary state, and then measured the order parameter at intervals larger than 
some estimated autocorrelation time. 

The results from three-dimensional simulations are in quantitative agreement with the 
mean field analysis, thus confirming that the upper critical dimension for this class of systems 
is 2. 

It is now interesting to see how this analysis is modified in the two-dimensional case, 
the physical one, when in the renormalization group sense only logarithmic corrections are 
expected. We simulated systems on a square lattice with linear size L up to 60. A cut of the 
phase space is shown in Fig. |7]. It is easily seen how the main features of the phase diagram 
remain unaltered in 2 dimensions; namely we can still identify a first order phase transition 
with a coexistence region where the system displays hysteresis (Fig. ||), and a crossover to 
a second order transition for large values of the noise. 

Similar results are obtained even in 1 dimension. It must be noticed how the coexistence 
region becomes smaller as one goes to low dimensionality, but is still present. 



VI. DISCUSSION AND CONCLUSIONS 

In a sense the combined effect of the attracting potential and of the noise is similar to 
that of a superposition of quenched an thermal noise (cfr. ||13||). Namely in both situations 
we have a force opposing the growth that depends on the height of the surface, and a noise 
term that randomly modulates in time and space such interaction. This is why we can use 
the quenched noise terminology and define e. g. the velocity exponent. Of course it is the 
attractive nature of the Morse potential that gives the transition its first order character, 
and we expect similar results for other short range attractive potential as well. 

We have shown how a simple modification of the nonlinear part of the KPZ equation 
can describe the dynamics of a surface that interacts with the substrate via a short range 
interaction. The mean field analysis has revealed that a first order unbinding transition 
takes place for a large range of the control parameters, and that in a whole region of the 
phase diagram the system's stationary position depends on the initial conditions (hysteresis). 
Furthermore a crossover to a second order transition takes place for large values of the control 
parameters, and the exponent controlling the vanishing of the order parameter depends 
continuously on the strength of the noise. 

The same situation is recovered in simulations above and below the upper critical di- 
mension d c = 2 
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FIGURES 



FIG. 1. Morse potential including the driving force —rh. 

FIG. 2. Graphical solution of the self consistence equation for (3 = 1, a = 0.4. 

FIG. 3. Hysteresis curve obtained from mean field (MF) and 2-d simulations. Here close and 
distant indicate the initial position of the surface with respect to the substrate. 

FIG. 4. The mean field phase diagram for fixed a = 1. The continuous line corresponding to 
r c = cr 2 /2 merges with the dotted line r C2 at the tricritical point (x). For larger (3 the transition 
is second order (dashed line). 

FIG. 5. The mean field phase diagram for (3 = 1. The continuous line corresponding to 
r c = a 2 /2 merges with the dotted line r C2 at the tricritical point (+). For larger a the transi- 
tion is second order (dashed line). 

FIG. 6. The relation between (3* and a* in a log-log plot 

FIG. 7. A cut of the phase diagram for the 2-dimensional system at fixed (3 = 1. Close and 
distant indicate the initial position of the surface with respect to the substrate. 
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